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A  solution     u  =  u(r,t)     of   the  differential  equation 
2 

is   to  be  found  for     r  >  r  ,   t  _>  0,   satisfying  the   initial 
condition 

u  =  0  foi'     t  =  0,  r  >  r^ 

and  the  boundary  condition 

u  =  u^  for   r  =  r,,  t  >  0. 

o  o'    — 

Our  nain   interest    is   to  find   out   at  what    time,    t,   at   a 

given   place,   r,   a   certain  percentage    of   the    final  temperature, 

u    ,    is'  reached.      Also  socie    Information   about   the   quantity 

du 
3r 

proportional   to  the  heat  flux  Is  desired. 

To  replace   equation    (1)   by  a  difference   equation, 

«e   proceed  as  follows: 

Let  u„  „  =  u(r^,t_)   and  replace   1)  by 
n,m  n     ni  " 


TE 


m 


*^n-fl,m  "  %,m  -  ^'n,ia  "  '^n-l.m 

Oi 

^  ^n                        ^  ^n-1 

^<^^n^-^n-l^ 

+      ^  ^n       2l       '^n-H,m  ~  "n-l,m 
^  ^n   *  ^n   •        ^  ^ n  -^    -^^n-1      ' 


Where 


^  'n  =  ^n+l-^n   *    ^  *ia  = 


n+1  "n   »    '-  "m        "m+l  "  *:;i     * 


-2- 
Now,  using  the  following  relationship 


y  '  ^^n-1  -~^ 


^"  =^n-l^l  "-p,    •'•  ^n-1  =T^^n 


n 


and  from  3) 

-^       ^^n-l--7~~r  T-+T  ^n  -  TT-H  ^^n 
Therefore, 

We  now  substitute  4)  and  5)  In   2)  and  obtain 

6)  n,m+l   "  %tn     _     %+i,in   "  ^n.ra         ^n,m  "  %-l,m 


^  r„  u    . ,        -  u     , 

.     ^    n  n+l,m  n-l.in 


Letting         . 


ex    ^_l     «'™  N(^ 

V)  — 1 5     =     6    ,   we   substitute    In     6)   and   obtain 

And   collecting  terms: 


9) 
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%,m+l  =  C  %-l,n   Ml   -  2C)   u^^^   +  C   u^^^^^ 
where  c  = -^i-^i-tJi 

L: 


The  results  obtained   lij  setting     C  =  77     are   ;  ot  as  re- 
liable   as  thoae   v.  1th  smaller     C: 

(T    =  ^  #       -      =5:  coofficlent     C  =  .3,    and   for 

a  finer  approximation 

^  =  ^  »        "     =  10s        coefficient     C  =   .275. 

Graphs  were  drawn  for  some  of  the  surfaces  at  a  constant 
distance  fron  the  center  of  the  hole.  They  indicate  that 
for  a  fixed  time  the  temperature  decreases  exponentially 
as  we  move  away  from  the  center  of  the  hole  and  that  for 
a  fixed  surface  the  temperature  Increases  very  slowly  to 
the    temperature   maintained   at  hole. 

The  small  difference   observed  when   comparing  the 
results  of  both  calculations,    indicates    that   the   Interval 
for        (^     =10     can   be   safely  employed,    if  v.e   do   not   ven- 
ture  too  far. 

A  rough  estimation  of   the   number  of   values  necessary 
to   indicate    the    state   of   the    surface   whose  distance    is 
20  times   the   bole   radius   is 

500  values   or  1500  multiplications   for    -   =  5 

1000  values   or  3000  multiplications   for    -  =  10. 
We,    therefore,   recommend   that   serjpus    consideration  be 
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glven  to   the  evaluation  of   the  following  formula"*'  v;hlch 
represents   the  exact  solution: 

10)     u  =  u     -^^     /'"-o.uS     Jq^^"^  ^o^^o^^   -  Jo(^o^)^o<^">   .  ^ 

o         w      l/        e  n n u 

Jo    <  V   *  ^  (^o") 

where  J       is  the   Besael  function  of  order  zero 

o 

Y       Is  the  Neumann  function  of  order  zero 

r^     Is   the  radius  of   the  hole 
o 

u       Is  the   temperature   of  hole 

The     nuaerlcal  integration  of  the   Integral  appear- 
ing In   (10)     is  straightforward,   if  one  makes  use   of  the 
following  approximate   formula,   valid  for  sriall   £  : 


11)       /    [integrand   in    (10)]du    ^log  —  [arctan   (a  log—     +  C)]^ 
o 

where      a   =|  =  ,6366,      c  =  |   (IT-  log  ^  )   =  -   .07380, 

I've  used    (11)   In  evaluating    (10)  for  the   particular  case 

2 


r 


-  2r      ,     ex  t  ~  4r''  and   obtained     u  =.66u 


•a-     U.S.    Cars  lav;   and  J.C.    Jaeger   -  Some   Two-Dliuensloncl 

Problems   in  Conduction  of   Ifeat  with  Circular   Symmetry   - 
Proceedings   of   London  Mathematical  Society  -  Series   2, 
volume   46,    1940,   pp.    361-338,    in   particular,    pp.    367-368. 

«■-»•  For     0  <    £  _<     u  ,  we  may  use   Simpson's  rule. 


-5- 


CASE  I 


^^^m  _  1       ^^n  _  1 
C  =  .3 


3 


^    ^1     ^2     ^3     ^4     ^5 


t^     100  0 

t,      100  30.0  0 

to     100  42,0  9.0     0 

tj     100  49.5  10.2  2.V  0 

100  54.7  22.1  5.94  .81  0 

100  58.5  27.0  9.25  2.11    .24   0 

100  61.5  31.1  12.4  3.69  .73   .07    0 

100  63.9  34.6  15.4  5.42  1.42   .25   .02    0 

100  63.9  37.6  18.2  7.21  2.27   .53   .08    .01  0 

100  67.6  40.3  20.7  9.03  3.23    .92    .19    .03  0   0 

100  69.1  42.6  23.1  10.8  4.28   1.39   .36    .07  .01  0  0 

100  70.4  44.7  25.3  12.5  5.37   1.95   .58    .14  .03  0  0 
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CASS   II 

o^/l  t^ 

/^-n 

_  1 

(^r„)^ 

4          . 

C   =   .275 

10 

t\r    r^    r^^    r^  v^  r^     Tg 


100  0 

100  27.5    0 

100  39. 9   7,56  0 

100  47,5  14.4  2.08  0 

100  52.8  20.1  4.90  .57  0 

100  56.8  24,9  7.89  1.60  ,16    0 

100  59.9  29.0  10.8  2.95  ,51   .04    0 

100  62.4  32.5  13.6  4.42  1.05   .16   .01   0 

100  64,5  35.5.  16.3  6.02  1,73    ,36   .05   0   0 

100  66.3  3B,2  18.8  7,67  2,53    ,65    .12   .01  0  0 

100  67.8  40.6  21.1  9.32  3.43   1.02   .24   . 04  0  0   0 

100  69.2  42.7  25.2  10.9  4.39  1.47   .40   .08  .01  0   0 
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